Rapid ramps across the BEC-BCS crossover: a novel route to measuring the 

superfluid gap 
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We investigate the response of superfluid Fermi gases to rapid changes of the three-dimensional 
s-wave scattering length a by solving the time-dependent Bogoliubov-de Gennes equations. In 
general the magnitude of the order parameter |A| performs oscillations, which are sometimes called 
the "Higgs" mode, with the angular frequency 2Agap/ft, where Agap is the gap in the spectrum 
of fermionic excitations. Firstly, we excite the oscillations with a linear ramp of 1/a and study the 
evolution of |A|. Secondly, we continously drive the system with a sinusoidal modulation of 1/a. 
In the first case, the oscillations in |A| damp according to a power law. In the second case, the 
continued driving causes revivals in the oscillations. In both cases, the excitation of the oscillations 
causes a reduction in the time-averaged value of |A|. We propose two experimental protocols, based 
around the two approaches, to measure the frequency and damping of the oscillations, and hence 
Agap. 



I. INTRODUCTION 

Nearly forty years ago Volkov and Kogan investigated 
the response of a Fermi superfluid to a small initial per- 
turbation of the order parameter A They predicted 
weakly damped oscillations in |A| with the angular fre- 
quency 2Agap/S., where Agap is the gap in the spectrum 
of fermionic excitations. These oscillations originate from 
the threshold for the creation of fermionic excitations by 
pair-breaking following a time-dependent perturbation. 
The threshold causes a branch-type singularity in the re- 
lated response functions at the frequency 2Agap/fi, and 
consequently oscillations in |A| with the same frequency. 
Hence a detection of the frequency would give a direct 
measurement of Agap. There is no exponential damping 
because the singularity lies on the real axis. However, 
since the singularity is not a pole, but a branch-point, 
the amplitude of the oscillations decreases with time ac- 
cording to a power law. Volkov and Kogan found that, 
in the Bardeen-Cooper-Schrieffer (BCS) regime, where 
the fermionic excitations are created with hnite momenta 
p near the Fermi momentum, the amplitude decays as 
i" 1 / 2 [1]. More recently, Gurarie calculated that the am- 
plitude decays as i~ 3 / 2 in the Bose-Einstein condensate 
(BEC) regime, where the excitations are created near 
the point p = ||. Persistent oscillations have been 
found for very large initial perturbations and particular 
nonequilibrium initial states • Due to a formal anal- 
ogy between the gap phenomenon in superfluid Fermi 
liquids and mass creation by the Higgs mechanism in the 
theory of elementary particles, these oscillations are also 
known as the "Higgs mode" [f|-Q- Despite some inves- 
tigation into the response of a superfluid Fermi gas to 
a modulation of the scattering length HITS, the Higgs 
mode has never been observed in Fermi gases. Anal- 
ogous phenomena have been observed in superconduc- 
tors fill. \1% and for Bosons near the Mott-insulator tran- 
sition 1131 . 



It is the purpose of this paper to provide a stepping- 
stone between the past theory described above and future 
experiments. By solving the time-dependent Bogoliubov- 
de Gennes equations [Tj-[l8[ we investigate the appear- 
ance of the order parameter oscillations in superfluid 
Fermi gases, following rapid ramps across the BEC-BCS 
crossover. In the spirit of Volkov and Kogan [1], we be- 
gin with relatively slow ramps of the three-dimensional 
s-wave scattering length a from a = do to a\ from t = 
to t = t± m H/Ef. We confirm that the angular fre- 
quency of the oscillations is given by 2Agap/ft, in which 

Agap = | A | is the BCS regime and Agap = \] + A 2 in 
the BEC regime, where /i is the chemical potential. The 
oscillations damp as predicted in Refs. [I], For abrupt 
ramps such that t\ <C h/Ef, we find that the order pa- 
rameter oscillates around a new value Aoo, which is less 
than the equilibrium value for a = a%. Note that our 
abrupt ramps never take the order parameter near zero, 
as instead studied in Ref . [5j . The frequency of the oscil- 
lations is set by the final value of the gap, which in turn 
depends on Aqq. The value of A^ decreases as we in- 
crease \l/a — l/ai\, see also Refs. [l9,l2d|- Furthermore, 
Aqo depends only on the instantaneous value of |A| just 
before the ramp, irrespective of whether the system was 
in equilibrium or not. We also study the response of |A| 
to a sinusoidal modulation of 1/a. In this case, |A| per- 
forms oscillations about a mean value which gradually 
decreases with t to a constant (see also Ref. |20|)- The 
oscillations do not damp, but show revivals due to the 
continued driving of the system. 

Building on these findings, we propose two experimen- 
tal protocols to detect the frequency and damping of the 
oscillations in |A| across the BEC-BCS crossover. The 
first is based on abrupt ramps of a, and the second on 
sinusoidal modulation of a. We confirm that the trap- 
ping potential in a real experiment would not signifi- 
cantly spoil the effectiveness of either protocol because 
the density of the gas does not respond on the timescale 
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of the oscillations in |A|. 

The paper is organised as follows. In section|TT]we out- 
line the system and methodology. In sections IIII Al and 
IIIIBI we investigate the response of a uniform Fermi su- 
perfluid to slow and abrupt ramps, such that t\ w h/Ef 
and t\ <C h/Ef, respectively. In section Ull CI we inves- 
tigate the response of a uniform Fermi superfluid to a 
sinusoidal modulation of 1/a. In section IlIIDI (jlll D 2|) 
we propose experimental protocol A (B), which is based 
on abrupt ramps of a (sinusoidal modulation of 1/a) . In 
section Qv] we conclude. 



II. SYSTEM AND METHODOLOGY 

We consider a three-dimensional superfluid Fermi gas 
with equal populations of the two spin components. We 
model its dynamics across the BEC-BCS crossover by 
solvin g th e time-dependent Bogoliubov-de Gennes equa- 
tions |l#jla |. Although this is an approximate theory, 
it is qualitatively correct and able to capture the main 
physics of the problem |21j. The equations are 
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where H = — fi 2 V 2 /2m + U(x) — [i, in which m is the 
atomic mass, U is the external potential and /i is the 
chemical potential. The order parameter is calculated as 
MM) = -gJ2n u r) v v> m which 9 is given by 1/fc/a = 
8TrE f /(gk 3 f ) + y/4EQ \n 2 E f ) [p. Here E f = h 2 kj/2m 

1/3 

and kf — (37r 2 n) are the Fermi energy and momen- 
tum of an ideal Fermi gas of density n, respectively. The 
cut-off energy E c is introduced in order to remove the ul- 
traviolet divergences in the Bogoliubov-de Gennes equa- 
tions with contact potentials. The density of the gas is 
n(r,t)=2J2 n \v ri (r,t)\ 2 . 

For most of the results in this article, we consider a 
uniform superfluid in a box with periodic boundary con- 
ditions and dimensions L x — 37.6kJ 1 and L± = 12.5kJ 1 
in the x- and other directions respectively. The poten- 
tial (7 = 0. Towards the end of the article, in order to 
simulate a more realistic system, we also consider the 
harmonic trapping potential U = muj 2 x 2 /2, where to x 
is the angular trapping frequency in the x— direction. 
Since, throughout the article, the potential U has no y 
or z dependence, we may write the functions u v (r, t) and 
v v (r,t) as u v (x,t)e t( - k y y+kzZ " 1 and v v (x, t)e^ kyV+kzZS > re- 
spectively, in which k y and k z are quantized according 
to k y — 27ta y /L± and k z = 2ita z /L±, where a y and 
a z are integers. We begin our simulations with station- 
ary solutions of Eq. (JT]) [l|| |22| for a — a . Figure Q] 
shows the equilibrium magnitude of the order parameter 
Aeq (1/fc/a) for a uniform Fermi gas. 
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FIG. 1: The equilibrium magnitude of the order parameter 
Aeq (1/kfo) obtained by solving the stationary Bogoliubov- 
de Gennes equations for a uniform Fermi gas. 



III. RESULTS AND DISCUSSION 
A. Slow ramps of the scattering length 

Our initial goal is to confirm the predictions of Volkov 
and Kogan [l| and Gurarie [2| for the period and damp- 
ing of the oscillations in |A|. In the spirit of Volkov 
and Kogan, we study the response of the superfluid to a 
small linear ramp of \/kfa from 1/fc/ao to 1/fc/ai over 
the timescale t\ « h/Ef. Both the small change in 1/a 
and the relatively large t\ minimise the reduction in |A|, 
specifically Aeq — Aoo, caused by the excitation of the os- 
cillation. In this work, we refer to these ramps as "slow" 
because they are close to being adiabatic. However, even 
these ramps are fast enough to be a challenge to realize 
in experiment. 

Figure E^a) shows the evolution of the magnitude of 
the order parameter |A(t)| in the uniform superfluid fol- 
lowing a ramp from l/kjao = —0.2 to 1/fc/di = (uni- 
tarity) over the time t\ — 2.2h/Ef. After the ramp |A| 
oscillates around the average value A^ = 0.664i?/ [indi- 
cated by the horizontal dashed line in Fig. Ufa)], which 
is very slightly lower (0.4%) than the equilibrium value 
at unitarity (see Fig. [1]): even this slow ramp of 1/fc/a 
causes some reduction in |A|. The period of the oscil- 
lation is A.7h/Ef. This agrees well with theprevious 
prediction of 7rft/Agap by Volkov and Kogan T], taking 
Agap = Aoo. We also quantify the damping of the os- 
cillation by evaluating the function A(t) = |A(t)| — A^ 
at the extremes of the oscillation, which are indicated 
by the crosses in Fig. [2ja). Assuming that A cx V , we 
plot log(.A) against log(i) in the Fig. [2ja) inset and mea- 
sure the gradient to determine 7. Our data is consistent 
with the prediction of Volkov and Kogan [l| , shown by 
the dotted line, that 7 and hence the gradient should be 
-0.5. 

We may gain a microscopic insight into the oscillations 
and their damping by studying the occupations of the 
energy levels O v = J \v v (r,t)\ dr. The occupations of 
the energy levels at the times indicated by the crosses in 
Fig. [21(a) are plotted as a function of |fc| with filled circles 
in Fig. [31 in chronological order from (a) to (f). Hence 
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FIG. 2: (a) The response of |A(t)| following a ramp from 
1/fc/Oo = —0.2 to l/k/ai — over the time ii = 2.2h/Ef. 
Inset: log(A) (see text), evaluated at the positions of the 
crosses in the main figure, against log(t). The dotted line is 
a fit with a gradient of —0.5. (b) As a, but l/fej>oo = 0.8, 
1/fc/oi = 1 and ti — 1.2h/Ef. The gradient of the dotted 
line in the inset is —1.5. The two insets are plotted with the 
same scale to aid comparison. 
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the left-hand (right-hand) column corresponds to maxi- 
mums (minimums) in |A|. (The occupations at the time 
indicated by the final cross in Fig.[2ta) are not plotted.) 
For comparison the dashed curve shows the equilibrium 
distribution function for 1/fc/a = 0. We see that the 
oscillations in |A| are due to a periodic flatening and 
steepening of the quasiparticlc distribution function. Af- 
ter an increasing number of oscillations, the distribution 
function becomes less smooth. This is the source of the 
damping of the oscillations. At large times this effect 
resembles a thermal excitation of the system, and hence 
< A e q (1/fc/a). 
To illustrate the corresponding behaviour in the BEC 
regime, Fig.rjjb) shows the response of |A(t)| to a ramp 
from 1/fc/ao = 0.8 to 1/fc/ai = 1.0 over the time t\ = 
1.2%/Ef. Following the ramp, |A| again oscillates around 
a new value A^ [indicated by the horizontal dashed line 
in Fig. H[b)], which is very close to but less than the 
equilibrium value for 1/k/a = 1 (see Fig.[TJ. We measure 
the period to be 1.6H/Et, which is slightly lower than 
the predicted 7rS/Agap = 1.7H/Ef, in which Agap is 
calculated as -y//x 2 + A 2 ^, taking /i to be the equilibrium 
value for 1/fc/a = 1. However, it is difficult to measure 
the period accurately because the damping is far more 
rapid than in Fig. [2ja). This is confirmed by the plot 
of log (.A) against log(t) in the Fig. [Z£b) inset. Our data 
points now follow the prediction of Gurarie Q , shown by 
the dotted line, that 7 and hence the gradient is —1.5 in 
the BEC regime. 



B. Abrupt ramps of the scattering length 

The slow ramps of the scattering length in the pre- 
vious section excite the oscillations in |A| with a min- 
imal reduction in the time-averaged |A|. However, the 



FIG. 3: (a) - (f): The filled circles indicate the occupation of 
the energy levels O v = J \v v (r,t)\ 2 dr at the points indicated 
by the crosses in Fig. Efa), in chronological order from the 
first maximum. Hence the left-hand (right-hand) column cor- 
responds to maximums (minimums) in jA|. The dashed curve 
shows the equilibrium distribution function for comparison. 



amplitude of resulting oscillations in |A| is quite small, 
and so it would difficult to base an experimental pro- 
posal around this scheme. In this section we discuss the 
response of the superfluid to an abrupt change in the 
scattering length. Here "abrupt" means short compared 
to H/Ef. This can excite oscillations in A with a large 
amplitude, and leads us to an experimental proposal to 
detect them. 

Firstly, we consider abrupt ramps to 1/fc/ai = (uni- 
tarity). Figure|4|a) shows the response of A(t) following 
an abrupt ramp from 1/fcjao = 0.2 (solid curve) and 1.0 
(dashed curve). Both ramps excite oscillations in |A|, but 
the time-averaged value of the order parameter, Aqo, is 
smaller for 1/kfao = TO. In general, we find that Aoo is 
always less than the equilibrium value Aeq(l/fc/ai) (see 
Fig-HJ), and that Aoo is smaller for large ramps. Notice 
that the period of the oscillations in |A| is set by A^, 
not by Aeq. Consequently, the period of the oscillation is 
longer following the ramp from 1/fc^ao = 1.0. Also note 
that the damping is more "noisy" following an abrupt 
ramp, so it would be more difficult to measure 7 by this 
method. 

Now we consider abrupt ramps to l/kfa% — 1 (the 
BEC regime). Figure H^b) shows the response of |A(t)| 
following an abrupt ramp from l/kja = —0.5 (solid 
curve), (dashed curve) and 0.5 (dash-dotted curve). 
Again we see that a large ramp leads to a small A m , and 
that as the ramp tends to zero approaches Aeq (1). 
As before, the value of A^ sets the period of the oscil- 



4 



(a) 



(b) 




5 10 15 

tE f /h 



2 4 6 8 



FIG. 4: (a) Solid (dashed) curve: the response of |A(£)| 
to an abrupt ramp at t = from l/k/ao = 0.2 (1.0) to 
1/fc/ai = 0. (b) Solid (dashed/dashed-dotted) curve: as 
(a), but 1/fc/ao = -0.5 (0/0.5) and 1/fc/ai = 1.0. Dotted 
curve: response of |A(t)| to an abrupt ramp at t = from 
1/fc/ao = to 1/fc/ai = —0.5 and a second at t = 0.90h/Ef 
tO l/fc/02 = 1.0. 



lations. Since we are in the BEC regime (1/kfd = 1) 
at the end of the simulation, variations in A co /Aeq(l) 
cause variations in the condensate fraction, which may 
be detected through analysis of the bimodal distribution 
following a free expansion (23|. 

The dotted curve shows the response of |A(t)| to a 
sequence of two ramps: the first from 1/fc/ao = to 
1/fc/ai = —0.5 at t = 0, and the second to l/kfa% = 1.0 
at t = 0.9QH/ Ef. We chose this time so that instanta- 
neous value of |A| just before the second ramp is equal 
to A e q (—0.5). We find that Aoo is identical to that for 
the single abrupt ramp from the equilibrium superfluid 
at 1/fc/ao = —0.5 to 1/fc/ai = 1.0 (solid curve). This 
illustrates that the superfluid has no memory of its pre- 
vious dynamics before the abrupt ramp, and that A^ is 
a true measure of the instantaneous value of |A| just be- 
fore the ramp. Although we show just one example here, 
we have confirmed that this is true for other parame- 
ters. This is an important requirement for the validity of 
experimental protocol A, described in section UlIDl 



C. Sinusoidal modulation of the scattering length 

An alternative method to excite oscillations in |A| is 
to modulate 1/fc/a at the resonant angular frequency 
2Agap/^i. This is potentially easier to realize experimen- 
tally because 1/fc/a would change smoothly. In principle 
the gap Agap(i) is a function of time, but we always mod- 
ulate 1/a at a constant angular frequency w m , according 
to 



1/fc/a (t) = 1/fc/ao + (0.04/fc/)sin (oj m t) 



(2) 
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FIG. 5: The response of |A(t)| to a sinusoidal modulation of 
1/fc/a at the resonant angular frequency ui m — 2Agap (t = 0) 
[see Eq. ©] with 1/fc/ao = 1-0 (top curve), 1/fc/ao = 0.5 
(middle curve) and 1/fc/ao = 0.0 (bottom curve). 



in the BCS regime and Ag ap (0) = y A|q (l/fc/a ) + /i 2 
in the BEC regime. The amplitude of the modulation 
is 0.04/ kf <C 1/fc/ so that we probe the response of the 
superfluid within a narrow range of 1/fc/a. Moreover, a 
small modulation of 1/a is easier to realize experimen- 
tally. The driving continues until t = t m = YiQh/Ef. 

Figure 02a) shows the evolution of |A(i)| in response 
to this driving for l/fc/a = 1.0, 0.5 and 0.0. For 
1/fc/ao = 1-0 (top curve), the magnitude of the oscil- 
lations increases slightly during the first two oscillations, 
then plateaus. Further driving causes no further amplifi- 
cation of the oscillations. The time-averaged value of |A| 
decreases only very slightly (on the order of 1%). This oc- 
curs because a small variation of 1/a in the BEC regime 
(large 1/fc/ao), where the fermionic atoms are tightly 
bound in molecules, is not sufficient to break the pairs. 
The response of |A(i)| is very different for 1/fc/ao = 0.5 
(middle curve) and 1/fc/ao = 0.0 (lower curve). Now 
we see that the amplitude of the oscillations grows over 
about five oscillations (see also Ref. 20]), as the time- 
averaged value of |A| slowly falls. However, on longer 
timescales we see a beating of the oscillations, and the 
time-averaged value of |A| saturates at a constant value. 
The beating occurs due to an interference between the 
driving frequency and the frequency of the oscillations in 
|A| set by A gap . 



D. Experimental protocols 

1. Protocol A 



In this section we drive at the resonant frequency by set- 
ting uj m = 2A g ap(0), where A gap (0) = A eq (l/fc/a ) 



Building on the abrupt ramps studied in section UlIBl 
we now propose an experimental protocol, which we call 
protocol A, to measure the frequency and damping of 
the oscillations in |A|. We perform two abrupt ramps 
on a superfluid Fermi gas with 1/fc/ao = 0.5. The first 
at t = excites the oscillations in |A| by ramping to 
1/fc/ai = 0. The second at t = t% ramps to l/fc/a2 = 1 
(the BEC regime). The dotted curve in Fig. [BJa) shows 
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FIG. 6: Experimental protocol A to detect the oscillations 
in |A|. (a) Dotted curve: response of |A(t)| in a uniform 
superfluid to an abrupt ramp at t = from = 0.5 to 

1/fe/Oi = and a second at t = t% = 7.1 h/Ef to l/fc/a2 
1 .0. Solid curve: as dotted curve, but for a 40 K superfluid in a 
trap with u x = 27r50 rad s" 1 and peak density 1.8 x 10 18 m~ 3 . 
Inset: profile of \A(x,t — 16h/Ef)\ in the trapped superfluid. 
The field-of-view is 120/t" 1 . (b) As (a), but with t 2 = 4.9 
h/Ej. (c) Aoo/Aeq in a uniform superfluid as a function of 
ta. 



the evolution of |A(i)| in a uniform superfluid for ti = 7.1 
h/Ef. In this case, ti coincides with a trough in |A|, and 
hence Aoo is small. The dotted curve in Fig. [(^b) shows 
the corresponding simulation for t 2 = 4.9 h/Ef. In this 
case, the second ramp coincides with a peak in |A|, and 
hence Aoo is large. 

In FiglHlJc) we plot A 00 /A e q(l) as a function of t 2 - 
The graph show damped oscillations that directly re- 
flect the behaviour of | A at unitarity following the initial 
ramp. As explained in section HUB I in the BEC regime 
Aoo/Aeq is a measure of the condensate fraction, which 
may be determined by ballistic expansion [23|, [24| . Hence 
this protocol will enable an experimentalist to resolve 
the frequency of the oscillations in |A| and their damp- 
ing across the BEC-BCS crossover. Of course, the ex- 
perimentalist would actually measure oscillations of |A| 
around a mean value which was slightly less than the 
equilibrium value. Hence the first ramp should be as 
small as possible, whilst maintaining a large enough vari- 
ation in Aqo to be measured in experiment. 

The solid curves in Figs. [BJa) and (b) show the results 
of the corresponding simulations for a K superfluid in 
a harmonic trap, with lo x = 2-7r50 rad s -1 and peak den- 
sity 1.8 x 10 18 m~ 3 . Although the oscillations in |A| 



are coupled to density oscillations in an inhomogeneous 
system, |A(f)| for the trapped superfluid is nearly iden- 
tical to that for the uniform superfluid because the den- 
sity cannot change significantly on a timescale of h/Ef. 
The insets show the profile of the order parameter at 
t = IQh/Ef. The profile in the Fig. UJa) inset has the 
smaller peak because t% coincides with a trough in |A|. 
Either side of the central peak the profile has some small 
shoulders, which occur because the frequency of the os- 
cillations in |A| decreases towards the edge of the cloud 
where the density and hence Agap is locally smaller. 
However, this is a minor effect which would only slightly 
reduce the variation in condensate fraction observed after 
expansion. 



2. Protocol B 

In this section we propose a second experimental pro- 
tocol, which we call protocol B, based on the periodic 
modulation of 1 /a studied in section IIII CI The advan- 
tage of this approach is that l/fc/a would vary smoothly, 
and hence protocol B would be easier to realize in exper- 
iment than protocol A. The disadvantage is that excita- 
tion of the superfluid saturates for large driving times t m , 
as explained in section Ull CI Hence it could be difficult 
to strongly excite the superfluid with this method and 
produce a large experimental signal. 

Protocol B is as follows. We excite a superfluid with 
an initial 1/fc/a = 1/fcjao using a periodic modulation 
of 1/a, as stated in Eq. |2]). At time t = t m we cease 
driving the superfluid and ramp to 1/kta — 1 over the 
time t m < t < t m 2- The timescale t m 2 — t m need not be 
short compared to h/Ef. The purpose of this ramp is to 
convert the excited superfluid with 1/fc/a = 1/fc^cto to 
an excited superfluid in the BEC regime, where the exci- 
tation may be quantified from the bimodal distribution. 
Hence, this final ramp may even be adiabatic. 

Figure [7] shows an example of protocol B with 
1/fc/ao = 0, t m = 107i7i/Agap(0) = A7.lh/E f and 
t m 2 —t m = 2.0h/Ef. The dotted curve in Fig.[7£a) shows 
the evolution of |A(t)| in a uniform superfluid for reso- 
nant driving [u) m — 2 Agap (0), see Eq. @], whilst the 
dotted curve Fig.[7fb) shows the corresponding result for 
off-resonant driving (oj m = 1.25 Agap (0)). Following the 
ramp into the BEC regime, the value of Aoo is much 
smaller for uj m = 2 Agap (0) , indicating that the reso- 
nant driving more strongly excites the superfluid. As 
explained previously, variations in Aoo/Aeq in the BEC 
regime may be detected by measuring the condensate 
fraction from the bimodal distribution following a free 
expansion. In Fig. EJc) we pl°t Aoo/Aeq (1) for a uni- 
form superfluid as a function of uj m . The minimum about 
LJ m = 2Agap (0) indicates the greater excitation of the 
gas at the resonant frequency of modulation. 

The solid curves in Figs. EJa) and (b) show the results 
of the corresponding simulations for a 40 K superfluid in a 
trap, with lu x — 2ir50 rad s _1 and peak density 1.8 x 10 18 
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FIG. 7: Experimental protocol B to detect the oscillations in 
|A|. (a) Dotted curve: response of |A(t)| in a uniform super- 
fluid to a resonant periodic modulation of a [cj m — 2Agap (0), 
see Eq. ((2])] followed by a ramp to 1/fc/a = 1 from t = 
t m = 47.1 to t = t m 2 = 49.1 h/Ef. Solid curve: as dotted 



curve, but for a K superfluid in a trap with ui x 



2tt50 



rad s _1 and peak density 1.8 x 10 18 m -3 . Inset: profile of 
\A(x,t = 60H/Ef)\ in the trapped superfluid. The field-of- 
view is 120k J 1 . (b) As (a), but for an off-resonant modula- 
tion of a [cj m = 1.25Agap (0)]. (c) A^/Alq as a function of 



m~ 3 . The solid and dotted curves in Fig. [7£a) diverge 
slightly for t > A0h/Ef because the density profile of the 
trapped superfluid begins to respond to the excitation of 
A|. This is illustrated by the Figs. [TJa) and (b) insets, 
which show the profiles of \A(x,t — 6QH/Ef)\. The pro- 
file in the Fig.lTJa) inset has a central dip. This occurs be- 
cause the gas is most strongly excited in the center where 
the frequency of the modulation is resonant. Towards 



the edges of the cloud the density is lower, and hence the 
local Ag a p is lower, meaning that the frequency of the 
modulation is larger than the local resonant frequency. 



IV. CONCLUSIONS 

Using Bogoliubov-de Gennes theory we have explored 
how the order parameter A in a superfluid Fermi gas re- 
sponds to rapid ramps across the BEC-BCS crossover. 
We have studied both linear ramps and sinusoidal mod- 
ulations of 1/a, where a is the three-dimensional s-wave 
scattering length. In general |A| performs oscillations 
with the angular frequency 2 Agap / % about a mean value 
Aqo, which is less than the equilibrium value of |A|. 
These oscillations are sometimes called the "Higgs mode" 
due to a formal analogy between the gap phenomenon in 
superfluid Fermi liquids and mass creation by the Higgs 
mechanism in the theory of elementary particles [a-ia]. 
The oscillations damp according to a power law follow- 
ing a linear ramp of 1/a. Continuous modulation of 1/a 
causes revivals of the oscillations. 

The main purpose of this work is to propose experi- 
mental protocols to detect the oscillations in | A| in Fermi 
supcrfluids, and hence measure Agap. We have proposed 
two protocols, refered to as protocol A and protocol B, 
which are based on abrupt ramps of a and modulation of 
1/a, respectively. Protocol A requires the experimental- 
ist to change the scattering length on a timescale faster 
than h/Ef. This is challenging, but may be achieved 
by using coils with low inductance, or by minimising Ef 
with dilute clouds or heavy atoms such a Yb [25j . Pro- 
tocol B is potentially easier to realize in experiment, be- 
cause 1/kfd would vary smoothly. Protocol A has the 
advantage that it could quantify the damping of the os- 
cillations. However, it has the disadvantage that it would 
measure the Agap following an abrupt ramp, and con- 
sequently the measured Agap may be significantly lower 
than the equilibrium value. In contrast, Protocol B ex- 
cites the equilibrium superfluid with a small-amplitude 
modulation of 1/fe/a, and hence may accurately measure 
the Agap of the equilibrium gas. 
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